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Calculation  of  Convective  Heat  Exchange 
During  the  Thawing  of  Water-Saturated  Soil 


N.  V.  Ukhova 

All-Union  Correspondence 
Engineering  and  Construction  Institute 


Usually,  when  calculating  the  depth  of  thawing  of  water- 
saturated  soil  as  a function  of  time,  the  concepts  of  the 
theory  of  thermal  conductivity  are  used.  Convective  heat  ex- 
change is  not  taken  into  account.  However,  in  large-skeletal 
soils,  and  particularly  in  rock  fill,  the  influence  of  con- 
vective heat  exchange  on  the  thawing  rate  may  be  considerable. 

V.  G.  Gol'dtman  (1958),  G.  V.  Porkhayev,  and  R.  M.  Sarkisyan 
(1960)  attempted  to  take  convective  heat  transfer  into  account 
in  studying  the  overall  heat  exchange  process  in  thawing  soil. 
They  obtained  an  approximate  solution  to  the  problem  of  heat 
transfer  caused  by  filtration  of  water  under  the  influence 
of  the  head  gradient,  arising  from  the  difference  in 
water  density  at  different  temperatures.  They  also  provided 
a joint  solution  for  the  equations  of  thermal  conductivity 
and  convection.  The  goal  of  the  present  paper  is  a develop- 
ment of  these  solutions  and  a comparison  of  the  results  ob- 
tained with  a different  statement  of  the  problem,  taking  into 
account  convective  heat  exchange  during  thawing  of  water- 
saturated  soils. 

Many  researchers  feel  that  the  movement  of  water  in  thinly 
dispersed  soils  obeys  Darcy's  Law,  while  in  coarse,  sandy 
gravel  and  other  materials  and  especially  in  fill  made  of 
boulders  and  rocks,  Darcy's  Law  is  not  applicable  [Aravin, 
Numerov  (1953);  Isbash  (1931)].  Since  this  problem  is  of 
critical  importance  for  the  processes  under  consideration, 
we  should  like  to  quote  the  opinion  on  this  subject  expressed 
by  Corresponding  Member  of  the  USSR  Academy  of  Sciences  M. 

A.  Velikanov  (1945);  "Beginning  with  the  analysis  of  the 
problem  as  a whole,  we  consider  it  particularly  necessary  to 
point  out  one  basic  defect  in  its  statement  in  the  literature. 
It  is  agreed  that  the  criterion  for  applicability  of  Darcy's 
Law  shall  be  the  dimensions  of  the  particles,  but  this  is 
not  absolutely  correct.  A linear  law  of  resistance  (and 
Darcy's  Law  is  one  of  these)  is  valid  according  to  the  gen- 
eral conditions  of  similarity,  applicable  in  modern  hydro- 
dynamics, at  small  values  of  the  Reynold's  number.  Conse- 
quently, and  with  respect  to  Darcy's  Law,  we  can  speak  of 
its  applicability  not  in  the  case  of  finely  particulate  soil, 
but  for  small  values  of  the  Reynold's  number."  And  also: 


>1 

< 


I 

( 


1 


"With  large  heads,  even  in  finely  granular  soils,  the  Reynold's 
number  may  turn  out  to  be  too  large  to  allow  Darcy's  Law  to 
be  used;  vice  versa,  in  coarse-grained  soils  with  very  small 
heads  (or  slopes),  such  small  velocities  may  be  produced  that 
Darcy’s  Law  turns  out  tc  be  quite  applicable." 


Figure  1.  Calculation  Diagram.  (At  top: 
thawed  soil;  at  bottom:  frozen  soil.) 


Since  a head  gradient  with  natural  convection  arises  from  the 
difference  in  the  density  of  water  at  different  temperatures 
and  is  consequently  very  small,  we  will  assume  that  convective 
mass  transfer  is  described  by  the  process  of  laminar  filtra- 
tion. 

By  analogy,  as  pointed  out  by  G.  V.  Porkhayev  and  R.  M.  Sar- 
kisyan, let  us  examine  the  following  calculation  diagram 
(Figure  1).  In  the  thawing  layer  of  soil,  at  some  depth  h 
from  the  surface  there  is  an  isotherm  t2  * 4°C.  Below  it,  at 
a depth  x,  is  the  boundary  between  the  frozen  soil  and  the 
temperature  ti  » 0°C.  From  surface  AA  to  surface  BB,  heat 
transfer  takes  place  as  a result  of  thermal  conductivity  as 
well  as  convection.  The  water  mixes  along  the  trajectory 
shown  in  Figure  1:  the  movement  of  the  water  downward  from 
AA  to  BB  unavoidably  causes  an  equal  volume  to  rise  from  BB 
to  AA.  The  soil  thawing  process  takes  place  quite  slowly,  so 
that  it  may  be  assumed  to  be  quasistationary. 

Looking  now  at  the  hydrostatic  pressure  at  levels  AA  and  BB,  we 
can  show  that  the  specific  motive  force  causing  mass  transfer. 
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as  a result  of  convection  arising  from  the  difference  in  the 
specific  layers  of  water,  will  be: 


(1) 


The  force  p,  causing  the  movement  of  the  water,  is  actually 
the  loss  of  pressure  along  the  length  of  the  filtration  path- 
way. Then  the  head  losses  in  this  section  will  be  expressed 
as  follows: 


where 


rip  J-ytp 


Y LtJj-. 

2._  __ 


(2) 


The  hydraulic  gradient  which  relates  to  the  length  of  the 
filtration  path  2x  will  be: 


LH  _ Yj,— y, 
~ 


Looking  at  filtration  in  an  ideal  soil,  where  all  of  the 
"current  tubes"  are  assumed  to  be  cylindrical  and  mutually 
parallel,  and  applying  Darcy's  Law,  we  have* 


itit 


Kt 


(4) 


where  wi  is  the  average  filtration  rate, 

K is  the  filtration  coefficient  of  the  soil,  and 
p is  the  porosity  of  the  soil. 


By  comparing  (4)  and  (3)  we  have: 


(5) 


* The  usual  form  of  Darcy's  Law  reads:  wt  * KI  (Editor's 
note . ) 
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Formula  (5)  is  valid  for  stable  filtration  when  the  water  is 
flowing  freely  through  a cross  section  BB.  At  the  same  time, 
end  BB  is  impermeable  to  water,  and  the  water  will  move  up- 
ward as  shown  above.  In  this  case,  due  to  the  movement  of  the 
water  in  the  opposite  direction,  the  average  filtration  rate 
will  decrease  relative  to  that  determined  by  (5).  In  order 
to  calculate  the  reverse  movement  of  the  water,  it  is  neces- 
sary to  introduce  a reduction  coefficient  n < 1 into  (5) . 


The  maximum  value  for  coefficient  n is  established  from  the 
following  conditions.  Half  of  all  of  the  "current  tubes"  in 
the  soil  will  be  occupied  by  water  moving  from  AA  to  BB  and 
half  by  water  moving  from  BB  to  AA.  This  will  halve  the 
actual  rate  of  movement  of  the  water  in  the  soil  relative  to 
the  average  filtration  rate.  In  addition,  due  to  the  fric- 
tion against  the  walls  of  the  tubes , the  water-velocity  dis- 
tribution over  the  cross  section  of  the  tube  will  be  in  the 
form  of  a paraboloid  (Figure  2) . In  this  case,  due  to  the 
formation  of  stagnant  zones  along  the  periphery  of  the  tube, 
the  average  velocity  of  the  water  through  the  entire  cross 
section  will  be  halved  again.  Herce,  the  maximum  value  of 
coefficient  n will  be  assumed  to  be  0.25  and  (5)  will  have 
the  following  form 
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Figure  2.  Diagram  of  Movement  of  Hater  in  Current  Tubes. 
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If  we  know  the  average  filtration  rate  of  the  water,  we  can 
determine  the  convective  heat  flux: 


-=CpO,yc#,  (^2  £i) 


(7) 


where  Cp  is  the  specific  thermal  capacity  of  water. 
Substituting  (6)  in  (7),  we  will  have: 

9l  = 0fi6i5^Cp(u-1,  (S) 

joint  computation  of  heat  transfer  caused  by  convection 
hermal  conductivity,  we  can  use  the  concept  of  the  ef- 
•e  coefficient  of  thermal  conductivity,  X3,  assuming  that 
the  heat  transmitted  by  contact  (in  other  words,  con- 
ation and  thermal  conductivity)  is  transmitted  only  by 
thermal  conductivity  (Kirpichev,  Mikheyev,  Eigensohn,  1940). 
Then  will  be  expressed  as  follows: 

hl  ~ t2 -6  (9) 


where  q = qi  + q2  is  the  total  heat  flux  due  to  convection 
and  thermal  conductivity. 

The  coefficient  of  increase  in  thermal  conductivity  due  to 
the  presence  of  convective  flux  will  be: 


(10) 


The  heat  flux,  determined  by  the  thermal  conductivity,  will 
be: 


where  X is  the  coefficient  of  thermal  conductivity  of  the  soil 
in  the  thawed  zone. 

Taking  (8),  (9),  and  (11)  into  account,  (10)  may  be  written 
as  follows: 


(12) 


BB1  AVAILABLE  COPY 


Formula  (12)  makes  it  possible  to  consider  jointly  the  heat 
transfers  due  to  convection  and  thermal  conductivity. 


The  value  of  the  coefficient  $ may  also  be  determined  on  the 
basis  of  a completely  different  statement  of  the  problem  than 
that  discussed  above.  Assume,  as  in  Figure  1,  that  the  total 
heat  flux  is  moving  from  surface  AA  to  surface  BB  (in  other 
words,  taking  convection  and  thermal  conductivity  into  ac- 
count) with  an  intensity  (Aerov,  Umnik,  1951) : 


(13) 


where  G ■ w x y&v  is  the  weight  velocity  of  the  water. 


The  general  solution  of  (13) , taking  into  account  the  boun- 
dary conditions  shown  in  Figure  1,  will  be: 


Then,  in  accordance  with  (9)  and  (10),  we  can  write: 


(15) 


The  movino  force,  as  before,  is  determined  by  (1).  On  the 
other  hand,  with  a stationary  regime  for  the  movement  of  the 
flow,  this  force  must  be  equal  to  the  hydraulic  resistance  of 
the  granular  layer  on  path  2x.  With  a laminar  regime  of  move- 
ment, the  hydraulic  resistance  of  the  granular  layer  will  be 
expressed  as  in  M.  E.  Aerov's  and  N.  N.  Umnik's  paper  (1951): 


(16) 


where  g is  the  acceleration  due  to  gravity, 

a is  the  surface  of  the  grains  in  a unit  volume  of  the 
layer,  and 

f is  the  coefficient  of  friction  of  the  granular  layer. 


The  coefficient  of  friction  for  low  flow  velocities  (with 
Re  < 20)  is: 

, % 36  ' a ^ * 

: ll7! 

where  u is  the  dynamic  coefficient  of  viscosity  of  the  water. 

Equating  (1)  and  (16)  and  taking  (17)  into  account  as  well, 
we  will  have: 

1 ^ " 


The  surface  of  grains  per  unit  volume  of  the  layer  may  be 
expressed  by  the  hydraulic  radius  of  the  "current  tubes"  of 
granular  layer  r [8]: 

n9> 

In  turn,  the  hydraulic  radius  may  be  linked  with  the  coeffi- 
cient of  filtration  of  the  granular  layer  through  the  coef- 
ficient of  permeability  [9] : 


where  8 s 1.87  is  the  parameter  which  takes  into  account  the 
shape  of  the  "current  tube"  in  the  soil. 
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From  (15)  and 


(22)  we  will  have: 
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(23) 


Formula  (23),  like  (12),  makes  it  possible  to  calculate  the 
combined  heat  transfer  due  to  convection  and  thermal  conduc- 
tivity in  thawing  soils.  The  formulas  obtained  are  approxi- 
mate and  based  upon  various  assumptions.  This  apparently  also 
explains  the  slight  difference  in  their  structures.  However, 
the  results  obtained  by  (12)  and  (23)  are  similar. 

If  these  formulas  are  converted  to  the  form 


® — | -0,0625— A 

• J n 


(12'  ) 


! 

where 


© — I =0, /94  A 

A — ^~P^*  V’») 


(23'  ) 


it  will  be  easy  to  compare  the  results  obtained.  Table  1 
shows  the  values  for  the  coefficients  at  A accordng  to  (12’) 
and  (23')  as  a function  of  the  porosity  of  the  soil  n. 


TABLE  1 


Porosity 

Formula  " 

0.3 

0.4 

0.5 

0.6 

12' 

0.208 

0.156 

0.125 

0.104 

23* 

0.194 

It  seems  to  us  that  the  structure  of  (12')  reflects  real  heat 
exchange  conditions  better  in  the  soil,  since  in  this  case  the 
coefficient  of  the  increase  in  thermal  conductivity  due  to  the 
presence  of  convective  currents  also  depends  upon  the  coeffi- 
cient of  porosity  of  the  soil. 
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Hence,  if  we  know  the  coefficient  of  filtration  of  the  soil 
in  the  thawing  zone,  we  can  calculate  the  value  <p.  When  <j>  < 

< 1.1  to  1.2,  the  influence  of  convection  on  general  heat  ex- 
change will  be  insignificant  and  the  rate  of  thawing  of  the 
soil  will  be  determined  by  thermal  conductivity.  When  1.2  < 

< <J>  < 4 , the  influence  of  convection  becomes  considerable  and 
it  must  be  calculated,  taking  into  account  the  thawing  rate  of 
the  soil. 

Wnen  $ > 4 convection  predominates;  in  other  words,  the  thaw- 
ing rate  of  the  soil  is  determined  not  so  much  by  thermal 
conductivity  as  by  convection. 

When  it  is  necessary  to  calculate  convection,  it  is  important 
to  calculate  the  value  of  the  effective  coefficient  of  thermal 
conductivity;  the  rate  of  thawing  of  the  soil  must  be  calcu- 
lated using  thermal  conductivity  formulas,  substituting  the 
values  obtained  into  them. 

For  the  sake  of  illustration  to  show  the  influence  of  con- 
vection upon  general  heat  exchange,  we  present  below  calcu- 
lations for  three  types  of  soils,  with  the  characteristics 
shown  in  Table  2. 


Granulometric  Composition 

cp 

cal/g 

degree 

A 

No. 

Soil 

0.05 

0.05- 

1.00 

1.00- 

3.00 

3.00- 

20.0 

20-80 

h 

cm/sec 

cal/cm, 

sec. 

degree 

1 

Fine  sand 

- 

23.5 

26.5 

50 

- 

3.  384 

0.0091 

0.404 

0.00628 

2 

Medium  sand 

- 

0.75 

14.0 

83.25 

- 

3.393 

0.098 

0.413 

0.0058 

3 

Gravel  and 
sand 

1.5 

10, 

,5 

46 

42 

0.35 

0.406 

1 

0.408 

0.0060 

The  characteristics  of  soils  1 and  2 are  taken  from  a paper 
by  P.  A.  Bogoslovskiy  (1956).  Characteristic  3 of  the  soil 
is  taken  from  the  work  of  V.  G.  Gol'dtman  (1958).  The  values 
of  A,  Cp  and  n are  missing  from  Gol'dtman "s  paper,  and  are 
included  conditionally  on  the  basis  of  corresponding  charac- 
• teristics  of  similar  types  of  soil. 
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T^e  values  of  4>  - 1 calculated  by  (12)  were  as  follows: 

Soil  1 p - 1 = 0.0000123  x 
Soil  2 p - 1 = 0.000144  x 

Soil  3 p - 1 = 0.000642  x 

where  x is  the  distance  between  the  +4°  and  0°  isotherms  in 
centimeters. 

Convective  heat  transfer  during  the  thawing  of  the  soils  in 
question  begins  to  be  expressed  (in  other  words,  p > 1.2)  at 
a distance  x equal  respectively  to  the  following: 

Soil  1 x = 162  m 

Soil  2 x = 13.9  m 

Soil  3 x = 3.12  m 

Hence,  thawing  of  the  first  and  second  types  of  soil  will 
occur  practically  as  the  result  of  thermal  conductivity  alone. 
The  thawing  of  the  third  type  of  soil,  at  a depth  of  more 
than  3 meters,  requires  that  convective  heat  transfer  also  be 
taken  into  account. 
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